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Abstract 

We apply factorization and eikonal methods from gauge theories to scattering amplitudes 
in gravity. We hypothesize that these amplitudes factor into an IR-divergent soft function 
and an IR-finite hard function, with the former given by the expectation value of a product 
of gravitational Wilson line operators. Using this approach, we show that the IR-divergent 
part of the ra-graviton scattering amplitude is given by the exponential of the one-loop IR 
divergence, as originally discovered by Weinberg, with no additional subleading IR-divergent 
contributions in dimensional regularization. 



1 Research supported in part by the NSF under grant PHY-0756518 
2 Research supported in part by the DOE under grant DE-FG02-92ER40706 
naculichObowdoin. edu, schnitzrObrandeis . edu 



1 



1 Introduction 



Infrared divergences afflict loop-level scattering amplitudes in both gauge and gravity theories. 
These IR divergences (as well as any UV divergences) may be regulated using dimensional regu- 
larization in d = 4 — 2e dimensions. In renormalizable theories, the UV divergences are removed 
from the renormalized amplitudes. The remaining IR divergences, which result from the exchange 
of soft virtual gauge bosons or gravitons, cancel against IR divergences arising from the emission 
of an arbitrary number of soft gauge bosons [IH3] or gravitons [HE] in inclusive cross-sections. 

In the scattering amplitudes of electrons or quarks, the leading IR divergence at L loops 
is proportional to l/e L . In the scattering amplitudes of massless partons, such as gluons (or 
massless quarks), the IR singularities are more severe, including both soft and collinear (or 
mass) divergences. In this case, the leading IR divergence is proportional to l/e 2L . In graviton 
scattering amplitudes, however, potential collinear divergences cancel from amplitudes due to 
momentum conservation [4] so that the leading IR divergence remains proportional to l/e L [6] . 

IR divergences exponentiate in both gauge [7] and gravity [4] theories, so that the leading 
L-loop IR divergence is expressed in terms of the leading IR divergence at one loop. Sub leading 
IR divergences in gauge theories have a more complicated structure, but decades of study have 
produced a detailed understanding [8H21"] which continues to evolve [221430] . The IR divergences 
of color-singlet form factors involving massless quarks or gluons are completely characterized by 
the coefficients of the beta function together with two constants at each loop order, the cusp 
7^ L ) and collinear anomalous dimensions [17] . The IR divergences of scattering amplitudes 
involve in addition a soft anomalous dimension matrix r^ L ^ (in color space [IB]) at each loop 
order |2Q1[2T]. 

The proof of the UV-finiteness of N = 8 supergravity amplitudes through three loops [31] has 
renewed interest in the detailed structure of gravity scattering amplitudes. J\f = 8 supergravity 
has now been shown to be UV-ffnite through six loops, based on the assumption of an £7(7) 
symmetry [32], and the situation for L > 7 has been discussed, for example, in refs. [3"2Tf3"5] . 
Using unitarity techniques, scattering amplitudes in Af = 8 supergravity can be expressed in 
terms of a basis of scalar integrals [31j[3llll0] . Using explicit Laurent expansions for these scalar 
integrals [4T1[42] the two-loop four-graviton scattering amplitude was computed in refs. |43|I4"4"]. 
That result revealed the absence of subleading, i.e., 0(1/ e), IR divergences at two loops in four- 
graviton amplitudes in = 8 supergravity. It was conjectured in ref. [43J that subleading IR 
divergences also vanish for n-graviton amplitudes (n > 4) at two loops, and possibly at higher 
loops as well. 

In fact, Weinberg's classic paper [4] showed that one-loop IR divergences of graviton ampli- 
tudes exponentiate to all loop orders, and argued that no additional IR divergences can arise from 
higher-loop diagrams because soft gravitons couple with a factor of the energy of the graviton. 
Weinberg regulated his amplitudes using UV and IR cutoffs rather than dimensional regulariza- 
tion (which had not yet been invented), but his arguments would imply the absence of subleading 
(in the e-expansion) IR divergences. 

In this paper, we extend the factorization approach which has been used so successfully in 
gauge theories [21] to gravity theories, confirming Weinberg's results, and setting the stage for 
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corrections to the eikonal limit. We assume that n-graviton scattering amplitudes can be factored 
into a product of IR-divergent and IR-fmite functions, and that the (helicity- independent) IR- 
divergent function can be expressed in terms of a set of eikonal loop diagrams encoded in the 
expectation value of a product of gravitational Wilson line operators. This approach is then used 
to show explicitly that potential subleading IR divergences indeed vanish. First, we reproduce 
the known one-loop IR divergence of the n-graviton amplitude in dimensional regularization 
[BJ. We then show that the two-loop IR divergence is entirely determined by the one- loop IR 
divergence, with no subleading IR divergences. In particular, various two-loop diagrams, which 
in gauge theory do contribute to the two-loop cusp anomalous dimension, have vanishing IR- 
divergent contributions in gravity. This result relies on a dimensional argument, which can 
be extended to higher loops. We then show (heuristically) that higher-loop diagrams, other 
than those contributing to the exponential of the one-loop IR divergence, will be IR finite, and 
will not result in subleading IR divergences in the dimensionally-regulated n-graviton scattering 
amplitude. In addition to confirming this perhaps to-have-been-expected result, our methods 
demonstrate yet another parallel between gauge and gravity amplitudes, which have been of 
much recent interest. 

This paper is organized as follows. In section 2, we express gravity scattering amplitudes in 
a factorized form, and in section 3, we hypothesize that the IR-divergent part can be written in 
terms of a product of gravitational Wilson lines, or equivalently, a set of eikonal loop diagrams. 
In sections 4 and 5, we compute the IR-divergent contributions of the one- and two-loop eikonal 
diagrams respectively. In section 6, we argue for the absence of subleading IR divergences at 
higher loops, and in section 7, we present our conclusions. 

Note added: After the first version of this paper appeared on the arXiv, White [15] applied 
path integral resummation techniques from gauge theories [16] to gravity amplitudes, finding 
strong supporting evidence for the factorization into soft and hard functions, and confirming 
the form of the soft function as a product of gravitational Wilson lines hypothesized below in 
eq. f)3.3p . He also extended this approach to massive scalar external particles, as well as to 
next-to-eikonal order. 

2 Infrared divergences of gravity amplitudes 

Gauge theory scattering amplitudes factorize into a product of jet, soft, and hard functions 
(see, e.g., ref. [21]). These separate functions describe, respectively, the contributions of gluons 
collinear to the external lines, soft gluons exchanged between external lines, and the short- 
distance (infrared-finite) scattering process. 

We will assume that, up to the loop order to which they are UV-finite, gravity scattering 
amplitudes obey a similar factorization, except that a jet function is not required since collinear 
divergences are absent in gravity [I]. We write the n-graviton scattering amplitude as 

A n = S n ■ H n (2-1) 

where S n is an IR-divergent factor describing the exchange of soft gravitons between the n 
external lines, and H n is IR-finite. Expanding each of the quantities in this equation in a loop 
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expansion in powers of k 2 = 32ttG, 



S n — 1 + ^ ^ 



(2.2) 



L=0 



L=l 



L=l 



we have 



(2.3) 



Exponentiation of leading gravitational IR divergences [4j suggests that we rewrite the soft func- 
tion as 



S n = exp 



L=l 



XL) 



i.e. 



<?(i) 



,(1) 



,(2) 



(2.4) 



(2.5) 



where the leading IR-divergence at L loops is given by [s^] L /L!, and the contain subleading 
IR divergences (if any). 

In Af = 8 supergravity, the all-loop-orders MHV four-graviton amplitude is proportional to 
the tree-level amplitude [4T)],I4"7] 



A 



M^Af\ 



H 



where M[ L ^ are helicity-independent functions of momenta, and M^" is finite as e — >■ 0. Com- 
bining eqs. (12.31) and (I2.6p . one has 



M { 4 Lf} Af\ 



(Lf) 



(2.6) 



M, 



(i) 



(i/) 



Mi 2) - 



M. 



(i) 



,(2) 



(2/) 



M. 



(i/) 



(2.7) 



Both m| 1} and m] 2) may be expressed as scalar integrals [ID] whose Laurent expansions in e are 
known explicitly |2JII2] . (The explicit expansions for the scalar integrals contributing to higher- 
loop amplitudes [3TJE1] are n °t yet known.) Using these, it was shown in refs. [HJHl] that the 
IR-divergent part of the two-loop four-point amplitude M± is given exactly by the IR-divergent 
part of one-half the square of the one- loop amplitude ; in other words, that sf^ is IR-finiteH 
In ref. [13], we conjectured that this result holds more generally for n-graviton amplitudes 
in M = 8 supergravity, i.e., that si 2 '' is IR-finite for any n. We also speculated that higher- loop 
subleading IR corrections, Sn for L > 2, might also be IR-finite. In the remainder of this paper, 
we present an eikonal calculation of the soft function S n which supports these conjectures. 



Explicit expressions for the IR-finite part of the two-loop amplitude were also obtained in refs. |43[l44]. 
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3 Eikonal description of the soft function 



In gauge theories, the soft function can be given an operator interpretation in terms of infinite 
Wilson lines (see, e.g., refs. [T9 | |2~1~ | |2"2"] ) . which can be expanded in terms of eikonal loop diagrams. 
Aybat, Dixon, and Sterman [21], employing results from ref. [48J, used an eikonal calculation to 
compute the two-loop soft anomalous dimension matrix for n-parton scattering amplitudes in 
QCD. We will adopt a similar approach to compute the soft function for gravity amplitudes. 

In ref. [J], Weinberg showed that the emission of a soft graviton (with momentum k) from 
an external line of an amplitude (with outgoing momentum p) multiplies the amplitude by the 
factoi0 

/ /,-\ Hi' i,' 

(3.1) 



pfj.pi> 



p ■ k + is 



(where k 2 = 32irG), independent of the spin of the external line. The numerator of eq. ( 13. ip de- 
rives from the vertex for graviton emission, and the denominator from the propagator 



(p + kf 



nr 



+ ie connecting the vertex to the rest of the amplitude, both in the eikonal 



limit k — > 0. 

We introduce a gravitational Wilson line operator 



$ p (a, 6) = Pexp [in ds p fl p u h fll/ (sp) 



(3.2) 



where s parametrizes a straight line segment x M = sp 11 . By Fourier-transforming the exponent 



of eq. ( 13. 2 p to momentum space, we see that it gives rise to a vertex in which a graviton 
of momentum k couples to a source (I3.ip . We hypothesize (in analogy with the gauge theory 
expressions in refs. [ZDI22]) that the soft function in the factorization ( 12 .ip may be expressed in 
terms of a product of infinite gravitational Wilson lines 



n$ Pi (o,oo; 



(3.3) 



The expansion of S n in powers of k then yields a series of eikonal loop diagrams in which 
soft gravitons are exchanged between n external lines of momenta pi. By our hypothesis, the IR 
divergences of these eikonal diagrams reproduce the IR divergences (both leading and subleading) 
of the original amplitude, as is the case for gauge theory. 

We observe that eq. ( 13. 3 p is the specialization of a more general Wilson loop operator intro- 
duced by Brandhuber, Heslop, Nasti, Spence, and Travaglini [H] (see also ref. [49J) 



W[C] 







Pexp I in (p dr x^(t)x v (r)h IJ , l/ (x(T)) 







(3.4) 



where x^(t) parametrizes a contour C. These authors explored whether MHV amplitudes in 
Af = 8 supergravity could be expressed in terms of gravitational Wilson loop expectation values. 



In this paper, we use the mostly-minus metric, ?yoo = 1, whereas Weinberg uses the mostly-plus metric. 
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(This was motivated by the discovery that the planar MHV n-point amplitude in M = 4 SYM 
theory, both IR-divergent and IR-finite parts, could be calculated in terms of expectation values 
of polygonal Wilson loops [50H56] .) As observed in ref. jH], however, the gravitational Wilson 
loop (13.41) . unlike the gauge theory Wilson loop, is not in general invariant under 

hp, ->• h^ + dptv + d v i p . (3.5) 

As was further observed in ref. |44j . if the contour C is composed of straight segments, the 
non-invariance of eq. (13.41) is concentrated at the endpoints of the segments. The cusps of the 
polygons result in gauge-variant expressions for polygonal gravitational Wilson loops, casting 
doubt on this prescription (although for a certain choice of gauge, the resulting expressions can 
be made to match the one- loop four-graviton scattering amplitude). 

In contrast, the variation of eq. (13. 3p under eq. (13. 5p vanishes because of momentum conser- 
vation 

n 

i=i 

provided that £ M (oo) is independent of direction. Consequently, the soft function S n is gauge- 
invariant. 

In the following two sections, we will use eq. (I3.3P to generate one- and two-loop eikonal 
diagrams to compute the IR divergences of n-graviton scattering amplitudes. 

4 One-loop IR divergences of gravity amplitudes 

Expanding the soft function (13. 3p to one-loop order gives rise to a set of eikonal diagrams, in 
which a graviton (thin line) is exchanged between an arbitrary pair i, j of external legs (thick 
lines), cf fig. [TJ 

We will find it convenient to evaluate the resulting eikonal diagrams in momentum space. 
Using the graviton propagator in de Donder gauge, 




Figure 1: The one-loop eikonal diagram contributing to S n ■ A graviton (thin line) is exchanged 
between external legs i and j (thick lines). The dots signify an arbitrary number of other external 
legs. 
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together with eikonal vertex (13. ip . we obtain 



where we have replaced k — > ft// 6 in eq. (13. 1 p (where \i is an arbitrary mass scale) to keep 
Sn dimensionless in d = 4 — 2e dimensions. We observe that the last term in the de Donder 
propagator (14.11) does not contribute to S n ^ (nor to any of the two-loop diagrams considered in 
the next section) because of the masslessness of the external lines, p? = p| = 0. 

The integral (14.21) formally vanishes due to cancellation between IR and UV divergences. The 
UV divergence is artificial (unrelated to any actual UV divergences of the gravity amplitudes), 
and results from the fact that we are using the eikonal (k — > 0) limit for soft-graviton exchange 
as if it were valid for all values of k. Once we remove these artificial UV divergences, Sn will 
give the IR divergence of the amplitude. 

We Feynman parametrize eq. (14. 2 j) as 

K^Vv^/n ^2 f dd k f°° , f°° , 1 



2 / U J ( 2n ) d J° Jo ' '[k* + 2k-(x lPi -x 2 p j 



is? 



K/i e Vr(l + e) , 2 f°° , f™ 1 



2 J ^Z^l dX2 [(x l p i -x 2 p j ) 2 ] 1+ ' ' ^ 

Since the external lines are massless (pf = p| = 0), this simplifies to 

W = (^) 2 ^^E(-^-^)" £ /^^/°°^(^)" 1 ~ £ - ( 4 - 4 ) 



i<3 



To disentangle UV and IR divergences, we follow the strategy of ref. [57] and use the symmetry 
under x\ x 2 to rewrite this as 



and then take e < to regulate the IR divergence at x 2 = 0: 



a?' = (f) 2 7S)^g(- 2 '"-») 1 " , (T)r' iXlirl " 2 '' (46) 

The integral over x\ formally vanishes as a result of cancellation of UV and IR poles. We cancel 
the UV pole by adding a counterterm — l/(2e) to the integral over x±, which leaves us with 

k^\ 2 r(i + e )^(-2 ft -p i ) 1 - 



n I 2 j (4vr) rf / 2 ^ 



e 



2 



2 



A x - (-2pi-pj) ( fi 
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where we have introduced 

\=^) 2 (4ne--Y . (4.8) 

The leading 1/e 2 pole in eq. (14.71) cancels via momentum conservation J2i<j Pi ' Pj = an d we 
are left with the IR-divergent contribution to the one-loop n-graviton amplitude 

This is precisely the result obtained by Dunbar and Norridge [6]. For n = 4, eq. (I4.9P reduces to 

+ 0(e°) (4.10) 

(where s, t, and u are the Mandelstam variables), which is indeed the IR-divergent part of the 
known j\f = 8 four-graviton amplitude [4T) | l4"3" t l4"4" l l58"]. 



S 



(i) 



s log 



+ tlog 



?) + " log V? 



5 Two-loop IR divergences of gravity amplitudes 

Expanding the soft function (13. 3p to two-loop order generates the same set of eikonal diagrams 
used in ref. [21] to compute the two-loop soft anomalous dimension matrix for gauge theories, cf 
fig. [2J It is again most convenient to evaluate these diagrams in momentum space. 

The ladder and crossed-ladder diagrams shown in figs. [2](a) and[2](b), respectively, give rise 
to the contributions 



K -f) E 



i<j 



E 

i<j 



d d h d d k 2 



(Pi ■ PjY 



(2n) d (2n) d kfk^h ■ Pj )(k 1 + fa) ■ Pj (h + k 2 ) ■ Pi [{h ■ Pi ) 



d d k 



-i{Pi - PjY 



(27r) d k 2 (k- Pi )(k- Pj ) 



(k 2 -Pi)_ 
(5.1 



The diagrams involving three eikonal legs shown in figs. Wis) an d E](d) contribute 



-\{ K 4) EEE 

j i^j k^i,j 



d% d d k 2 



{Pi-Pjf{Pj -Pkf 



(2vr) d (2yr) d k\k\{k x ■ Pl )(k 2 ■ Pk )(h x + k 2 ) ■ Pj 



*m EEE 

j i^j k^i,j 



d d h 
{2it) d 



-i(Pi - Pj] 



d d k 2 



_(k 2 - Pj ) 
-i(Pj -Pk) 



fciCfe -Pj)(h - Pk 



(h- Pj )_ 

(5.2 



The diagram involving four eikonal legs in fig. |S^e) contributes 



\m ee 



i<j , k < 1 . 



d d h 

(2n) d 



-i{Pi • Pj) 



d d k 2 

(2n) d 



-i(Pk - Pi) 2 
k\{k 2 - Pk ){k 2 - P i) 



(5.3) 





(h) 



(i) 



G) 



Figure 2: Two- loop eikonal diagrams contributing to Sn 



(2) 



The factors of | in eqs. (15. 2 p and H5 .31) compensate for double counting some of the diagrams. 
The sum of the Abelian-like contributions (I5.ip - fl5.3p can be neatly rewritten as one-half the 
square of the one- loop eikonal contribution (14. 2p . namely 



E 

Kj 



d d k —i(pi ■ Pj) 2 
(2vr) d k 2 {k -pi)(k ■ pj) 



(5.4) 



as was observed by Weinberg long ago [I]. 

The 1PI diagrams shown in figs. |2(f) and[Hg), in which one graviton is emitted and absorbed 
by the same eikonal line, contain factors of p 2 in their numerators, and hence vanish for massless 
external lines. 

Figure [2](h) involves a triple graviton vertex 

3 



E 

i<j<k 



d d h d d k 2 d% 



S^ih + h + h) 



(2n) d (27i) d (27r) d k\ k\ k\ (h ■ Pl )(k 2 ■ Pj )(k 3 ■ p k ) 



V 3 (k 1 ,k 2 ,k 3 ;p i ,p j ,p k ) (5.5) 



where 



V 3 (k l} k 2 , k 3 ;pi,pj,p k ) = PiPiPjPjPlplG^ap^ih, k 2 , k 3 ) 
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(5.6) 



is given by the numerators from the gravitational eikonal vertices contracted with the three- 
graviton vertex G tlUta ^ !(TT (ki, k 2 , k 3 ). This vertex is in general gauge- dependent. Using the ex- 
pression obtained by Bern and Grant! one has [59J 



V 3 (k 1 ,k 2 ,h;p i ,p j ,p k ) 



(Pi ■ Pj){Pk ■ h) + (Pj ■ Pk){Pi ■ k 2 ) + (p k ■ Pi){Pj ■ h) 



1 2 



+ 



(Pi ■ Pj){Pk ■ k 2 ) + (pj ■ Pk){Pi ■ h) + (p k ■ Pi){pj ■ h 



2 



(5.7) 



which is manifestly even under ki — > —ki. The three-graviton vertex in other gauges is more 
complicated (see e.g. ref. j5]) but shares the property of being even under ki — > —ki. This 
implies that the integrand in eq. (15.51) is odd under ki — > —ki, hence the integral vanishes 
independent of gauge. It therefore has no UV divergence, and no counterterm, so this diagram 
does not contribute to the soft function. (In nonabelian gauge theories, the analogous diagram 
also vanishes |21j.) 

Figure [2](i), which contributes 



i<j 



d d k x d d k 2 d d k 3 



5^(h + k 2 + k 3 



(27T) d (27T) d (27T) d k\ k\ k\ (h ■ Pi )(k 3 ■ Pl )(k 3 ■ p.j) 



V^k^k^k^piiPiiPj) (5i 



also vanishes because the integrand is odd under ki — > —ki. Hence, this diagram does not 
contribute to the soft function either. (In nonabelian gauge theories, the analogous diagram does 
contribute [2T].) 

The final diagram, shown in fig. Etj), contains a self-energy insertion in the graviton propa- 
gator exchanged between two eikonal legs, and contributes 



i<j 



d d k P^P^P^vMk) 
(2vr) d k ■ pi (k 2 ) 2 k ■ Pj 



(5.9) 



where the one-loop graviton polarization tensor has the form 



W v ^(k) 



k^^k^ik 2 ) + k 2 (k^k u r] af} + k a k^ u )U 2 (k 2 ] 

+k 2 (k^k a r] u(3 + k^r] va + k v k a rfP + k v k^rf a )Jl 3 {k 2 
+ (k 2 ) 2 7 1 ^7 1 a ' s U 4 (k 2 ) + {k 2 ) 2 {rf a r\ ufi + 7]^rj ua )Il 5 (k 2 ) 



(5.10) 



5 The Bern-Grant gauge was developed for calculating tree-level graviton scattering, and has not been validated 
for loop calculations. However, use of the Bern-Grant gauge does not alter any of the eikonal results computed 
in this paper provided the external lines are massless. 
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Contractions with the external momenta in eq. (15. 9p leaves 



~) h' « 



(p) 2 1V y p JV y (k- Pi )(k- Pj ) 

(5.11) 

Assuming the fields contributing to the graviton self-energy are massless, the Ui(k 2 ) all have 
momentum dependence (k 2 )~ e , so the integrands of each of the terms in eq. (15. lip go as (k 2 )^ 1 ^ 6 , 
and therefore the integral does not diverge in the IR region k — >■ 0. Consequently, fig. WS) can 
only give an IR-finite contribution to the n-graviton amplitude. 

To summarize, the five two-loop eikonal diagrams fig. EJa)-(e) add up to give precisely one- 
half the square of the one- loop eikonal diagram. None of the remaining five diagrams fig. Eff)-(j) 
contain IR-divergent contributions. Therefore we have 

Sf) = 1 [£«] 2 + (IR finite) ==* s {2) = (IR finite) (5.12) 

i.e., there are no subleading IR-divergences in the two-loop n-graviton scattering amplitude. This 
confirms the explicit finding for the two-loop four-graviton amplitude in M = 8 supergravity 
jHim], an d extends it to the n-graviton amplitude. 



6 Higher-loop IR divergences of gravity amplitudes 

In the previous section, we computed the eikonal contribution to the two-loop n-graviton scat- 
tering amplitude, and explicitly showed the absence of subleading IR divergences. These cal- 
culations, and in particular that of the self-energy diagram, suggest a strategy for showing the 
absence of subleading IR corrections at arbitrary loop order. While we make no claim to rigor, 
the arguments of this section suggest that the only diagrams contributing to the IR-divergent 
part of the amplitude are those generated by the expansion of the exponential of the one-loop 
diagram. 

Consider an L-loop n-graviton eikonal diagram containing E eikonal vertices at which a soft- 
graviton attaches to an external line. Schematically such a diagram gives a contribution 

n / d% f^ ava {£) n i- Kj f) f§ . (6.i) 

j=l J 0=1 ^ J Pa a 

Here, p a is the momentum of the external line from which the soft graviton is emitted; the 
p a are not necessarily distinct. The momentum of the emitted graviton, k a , is some linear 
combination of the independent loop momenta £j. The function f^ 1 u 1 --- fJ , E u E (£) is composed of 
internal propagators and vertices. Overall, the internal vertices contribute a factor of (kh € ) 2L ~ e 
in order that the whole diagram goes as (k//) 2L . Assuming that all the interacting fields are 
massless, f^ lUl -fM E u E (^) is homogeneous of degree — 2L in the loop momenta £ to make the final 
result dimensionless. Hence, the integral scales as 

M'VjiSs (6 2) 
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• • • 

• • •■ 

(a) (b) (c) (d) 

Figure 3: Some three-loop eikonal diagrams 

At this point, we divide the L-loop eikonal diagrams into two classes: those with E = 2L 
and those with E < 2L. The former class involves independent exchanges of L soft gravitons 
between external lines, e.g. figs. E]Ja)-(e) of the previous section, or fig. [3]^a). The integrals over 
loop momenta in eq. ( 16.21) generate IR divergences in these diagrams. By a generalization of the 
calculation of the previous section, however, the sum of all such Abelian-like diagrams is given 
by@] 

^m L . (6.3) 

(Diagrams in the class E = 2L in which a soft graviton is emitted and absorbed by the same 
external line, e.g. figs. [2]T) and (g), vanish before integration because of the presence of factors 
of p 2 in the integrand.) 

The diagrams in the class E < 2L all contain non-trivial internal vertices, e.g. figs. E^b)-(d). 
Diagrams like fig. EJ^b) also have IR divergences, but these are generated from the IR divergences 
of lower loop diagrams. For example fig. G](b) contributes to the term Sn^Sn^ in the expansion of 
S n in eq. (JZU). 

Finally, diagrams like fig. [3](c) and (d), in which all the soft gravitons are connected in a 
connected web [SUES], are not generated by lower- loop diagrams, and contribute directly to Sn . 
One may characterize such diagrams as non-Abelian-like. Included in this category are diagrams 
with vertices containing four (or more) gravitons. If the loop momenta in these diagrams are all 
scaled uniformly, power counting suggests that eq. (16.21) has no IR divergence as £ — > 0. This 
suggests that the functions sffl are all IR-finite. Consequently, subleading IR divergences are 
absent in the n-graviton amplitude. 

7 Conclusions 

In this paper, we have extended to gravity the powerful factorization methods used in gauge 
theory. These methods imply that scattering amplitudes can be factored into an IR-divergent 
soft function and an IR-finite hard function, and that the soft function can be written as the 
vacuum expectation value of a product of gravitational Wilson loops. Using this approach, we 
calculated the IR divergences of n-graviton scattering amplitudes, and confirmed the classic result 




12 



of Weinberg [4] that these are generated by the exponential of the one-loop IR divergences, with 
no additional IR-divergent contributions. 

There are various open issues which deserve mention: 

(1) In this paper we assumed that the factorization of amplitudes is valid, at least to the order 
in perturbation theory to which the gravity theory is UV finite. In ref. [IS], strong supporting 
evidence for this hypothesis has been presented. 

(2) What happens at the loop order at which the gravity theory is UV-divergent, and thus 
requires UV subtractions? Does the factorization fail, or do the regulated UV divergences and 
accompanying undetermined parameters become part of the definition of the hard function, with 
factorization intact? 

(3) What happens for massive external Wilson lines? White jU] has now studied this general- 
ization. 

It should also be mentioned that eikonal resummation in gravity and supergravity has been 
used in ref. [63] to study high-energy scattering in the ultra-Planckian limit. Although there does 
not appear to be a direct connection with our use of the eikonal approximation to determine IR 
divergences, this issue deserves further consideration. 
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